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ON THE EXTENSION OF DELAUNAY'S METHOD IN THE LUNAR 

THEORY TO THE GENERAL PROBLEM OF 

PLANETARY MOTION* 

BY 

G. W. HILL 

Part I. — Exposition of the Theory. 

The method of integrating the differential equations of motion, adopted by 
Delaunay for the elaboration of his lunar theory, as it is explained by him, 
demands its division into several cases, and is established through very tedious 
transformations. These disadvantages disappear when the greatest generality 
is given to the procedure. Hence, an explanation of the method, as it would be 
applied to the motion of a planetary system like the solar, will, doubtless, be 
welcome to astronomers. 

I. 

Let T denote the living force of the system, Q, the potential function, and, 
with Poincare, put 

F=D,- T. 

The k variables, necessary for completely defining the position of the system, 
may be denoted by ft , ft , • • • , q k . Use accents to denote complete differen- 
tiation with respect to the time of the latter, and we have 

r= function (ft, q 2 , •••, q k , ft , ft, •••, q' k ) . 

The partial derivatives of this function with respect to the k variables q\ are to 
be used as variables instead of the latter, and we put 

dT 
dq. 

By means of these k equations the q' { can be eliminated from T, and thus will 
result : 

T= function (ft , ft, •••, q k ,p,,p 2 , ■ ■ • , p J . 



* Presented to the Society February 24, 1900. Received for publication February 5, 1900. 
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Then the system of differential equations for determining the variables p. and 

q { is : 
* dp. BF do. BF 

K > dt dq^ dt Bp { ' K ' ' 

Let us suppose that D, is separated into the two parts O and O t , and that, 
when we neglect 0, t in F, the equations (1) can be completely integrated, and 
their integrals expressed in terms of two sets of h quantities each, symbolized 
thus : 

A > A > • * - ' A ' 
\i \i • • • i \i 

of which the first set are constants, and the second set linear functions of the 
time of the form n.t + c. , n. being a function of the L t , and c. an arbitrary 
constant. Nothing forbids our taking the L t such that they may be the ele- 
ments severally conjugate to the \ . 

Now, desiring to integrate the equations (1) when F has its complete value, 
we may adopt the A an( ^ *^ e \ as ^ e dependent variables to be employed. The 
differential equations of the problem are then : 

dL. BF dX. BF 

( 2 ) Tir-aT' -m-~nr.> c=i.».-.*)- 

Here the function F has been made to involve the L i and \ t by eliminating the 
old variables p. and q. from it by means of their values given by the integrals 
derived on the supposition that O = O . As 

F= a constant 

is an integral of the problem, and O — T= a constant, when il t is neglected, 
it is quite evident that when we substitute in il — T for p i and q. their values 
in terms of the A an( ^ \ > * ne \ completely disappear and fi — T becomes a 
function of the L i only. Thus, in the second form for F, the variables X ; 
enter into it solely through the portion O x . 

II. 

In order to exemplify we will adduce the solar system composed of the 
Sun and the eight major planets. We will suppose that the masses of the Sun, 
]\(Iercury, Venus, ■ • • , Neptune are denoted by m , m x , m 2 , • • • , m g , and will 

put 

m. 
Mi= m o + m i + m 2 + M ra i' K t = ^> (* = 0, 1, ■••, 8). 
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Let the type of representation of the rectangular coordinates of the planets rela- 
tive to the Sun be as follows : 

Mercury x x , 
Venus x 2 + k x x x , 
Earth x s + k 2 x 2 + k^ , 



Neptune x s + « 7 tc 7 + • • • + K t x t . 
The differential equations these variables satisfy are : 

d 2 x. da 



(3) 



'4—1 % 



df 



Ik 



dx. ' 

% 

da 



^-i*« df 

d\ 
^-^ df ~ 



da 



= 1, 2, 



8). 



dz. 



Here X2 denotes the sum of the products of every two masses of the system di- 
vided by their distance, a relation we will write thus : 



(4) 






mm. 

Jk 3 



Suppose that the portion to be separated from ft is 

(5) o>- «*!;?• 



r. standing for */{x\ + y'\ + »*) • Then, if ft is substituted for O in equations 
(3) , and the members are divided by M i _ 1 « i , we get 



(6) 



dt 



u. x. 



Jf 



y. 

i + m, 



a. , r: 

~i— l * 

' t t'i _____ A 

>. , r 3 ~ ' 



(<=1, 2, 



C?%. w. z. 



It will be seen that each group of these equations, corresponding to the same 
value of i , is independent of all the rest, and that it differs from the group of 
equations of relative motion of two bodies only in that the constant m /vX-i 
takes the place of m + m. . 



208 



G. W. HILL: EXTENSION OE DELAUNAY'S METHOD 



[April 



Let a. be the semi-axis major, e i the eccentricity, A. the inclination, I the 
mean anomaly, g. the angular distance of the perihelion from the node, and h, 
the longitude of the node. Put 



(?) 



\ ° a. , 8 ' 



IT,— G, cos A. , 



(* = 1, 2, ••-, 8). 



Then, when the elements become variable by reason of the addition of fl x to fl , 
they will satisfy the differential equations : 



f dL t dR. 


dl dR 


It ~"dJ. ' 


dt ~ ~ dL. 

1 


dG, dR. 
~dt = 'dg]' 


dg ( dR ( 
dt ~ dG. 


dH. dR. 


dh. dR. 


dt ~ dh. 


1t =: ~'dir. 



(8) 



where R i will be, in terms of O s , mentioned above, 

Pi 



1 = 1, 2, ■•-,8), 



(9) 



m„ 



2a. ^ /». .*. 



(* = l,2,---,8). 



Desiring to have the same perturbative function, whatever may be the integer i, 
we multiply the values (7) of L. , G t , H., as also the value (9) of R. by the 
constant f^^^, which does not alter the form of equation (8). We now have : 



(10) 



as also : 



H.= G^osA., 



(i=l,8,---,8), 



(") ^=».2a+».!>,[ s L -J]+rp. 

If the planetary coordinates in the last equation are replaced by the elements 
(10) and the I., g t , h., the differential equations of the system are : 
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(12) 



dL t 


dF 


dl. dF 


dt 


~ dl/' 


dt ~d~L. 

% 


dG ( 


dF 


dg { dF 


dt 


"¥,' 


dt ~ dQ. 


dH. 


dF 


dh t dF 


dt 


~ dh: 


di~~"dW. 



(i = l,2,---,8). 



In the second term of the right member of (11) the quantities 1/A 0)2 , 1/A 03 , 
• • • , 1/A^g, can be developed in infinite series, the first terms of which are 
l/r 2 , l/r s , • • • , l/r g , and thus are cancelled by the term 1/r. . Then the two 
latter terms of (11) are of the second order with reference to planetary masses. 

HI. 

In order to make the application of Delatjnay's method it appears necessary 
that F should be developed in a series, finite or infinite and periodic with respect 
to the variables L, g t , h i , which have been named the angular variables. 
In astronomical problems the series are generally infinite. For legitimate em- 
ployment this series must remain convergent throughout the whole duration of 
motion, while t is passing from — oo to + oo . It becomes then pertinent to ask 
what conditions must be fulfilled in order that this series may be convergent. 
It is well known that the reciprocal of the distance between two planets can be 
developed in a convergent infinite series, periodic with respect to the mean an- 
omalies of the planets, provided that the orbits, as they stand in space, have no 
point in common, or when the reciprocal of the distance never becomes infinite. 
The condition of convergence in the present case is precisely similar to this. 
Here, however, not only the mean anomalies I. are left indeterminate in the series, 
but also the remaining angular variables g. and h i which define the positions of 
the perihelia and nodes. Hence, in the present case, there must not only be no 
actual intersection of the orbits, but none when the perihelia and nodes are 
shifted in every possible way, the linear variables, or the mean distances, eccen- 
tricities and inclinations retaining their actual values. In the Delaunay develop- 
ment of the reciprocal of the distance between two planets, it is necessary and it 
suffices for convergence that the perihelion radius of one of the planets should 
always exceed the aphelion radius of the other. 

We may consider this subject under a more general aspect. Let F have the 
periodic development 

(13) F = ZA cos lj r \ + j 2 X 2 + ...+ j k \2 , 

where X l • • • \ are the angular variables, the j positive or negative integers, and 
J. is a function of the linear variables L only. That this infinite series may be 

Trans. Am. Math. Soc. 14 
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convergent, F must not only actually never become infinite, but never even 
potentially so. It is necessary here to explain what we mean by the qualifying 
epithet " potentially." If, while the linear variables L are supposed to maintain 
their actual values, and, consequently, the coefficients A their actual values, we 
allow to all the angular variables X the complete swing of movement from to 2w , 
F remains always finite, we say it never -potentially becomes infinite. In order 
that F may not actually become infinite it is necessary and sufficient that the 
velocities of and the distances between the points of the system should remain 
finite. In order that F may not potentially become infinite, it is necessary and 
sufficient that the values of the linear variables L x , L 2 , • • • , L k should remain 
within a certain domain. The definition of this domain is very complex after a 
Delaunay transformation has been operated, but is quite simple in terms of the 
original Keplerian linear variables L x , L 2 , • • • , L k . 

We may illustrate this subject by bringing forward the case of the solar sys- 
tem as it has been described in § II. Employing the linear elements a { and e. 
of equation (10) or of equations (7), the inequalities which define this domain are : 



(14) 



«! (1 - e,) > , 

°*(i - e 2 ) - (1 - *i)Oi(l + e i) > ° » 

«3 (1 ~ «j) ~ (1 ~ O % (1 + O > , 

«8( 1 - e 8)-( 1 -«r) a 8( 1 + e 7)> ' 
K(l+e 8 ) <oo, 



or, if we prefer to use the elements £/ i and G. of the same equations, the in- 
equalities are : 



(15) 



ADA - ^A -GU>o, 
£§ L,[L t - v'ZTTGa > £| ACA + ^A^?] , 

|§ i 3 [A - ^Li^GU > £f A[A + ^A^] . 



A[A - ^A - GS > 5 ADA + ^r - #?] . 



A+^-<72<. 



It is remarkable that the inclinations of the orbits play no part in these in- 
equalities. 



1900] TO THE GENERAL PROBLEM OF PLANETARY MOTION 211 

IV. 

As, in general, it is not necessary to distinguish between the three kinds 
I. , (j. , A i of angular variables, nor between the three kinds L i , G i , H { , of 
linear variables, for simplicity of notation we shall suppose that the angular 
variables are denoted by l x , l 2 , • • • , l k , and their corresponding conjugate 
linear variables by L l , L 2 , • • • , Z/ 4 . 

Selecting a particular linear combination of of the angular variables so that 

the J being positive or negative integers prime to each other, Delaunay's 
method, somewhat generalized, consists in making such a transformation of 
variables as would constitute a complete solution of the problem if F, in its 
periodic development, contained as arguments only integral multiples of 8 . That 
is, in this special case, the new linear variables would turn out constants, and 
the new angular variables would be of the form n(t + c) , n and c being likewise 
constants. It is clear that when we make such a transformation in F, the 
terms in the former periodic development involving the cosines of the finite 
multiples of will disappear, but the absolute term will receive a modification. 
A little consideration will make it evident that the derivation of such a trans- 
formation is dependent on quadratures. 

The discussion of this derivation is greatly facilitated by making a linear 
transformation of variables, in the cases of both angular and linear variables. 
In this it is evident that we can take 6 as one of the angular variables ; then let 
% be its conjugate linear variable. Thus, we may have the following as the 
variables involved in the problem : 

Linear variables, © , A t , A 2 , • • • , A t _j ; 
Angular variables, 6 , \ , \ 2 , • • • , X,.^ . 

And the canonical system of differential equations will be : 



(16) 



d® dF 


dd dF 


dt ~ ~d6 ' 


di~~ m ' 


a\ x dF 


d\ dF 


dt ~ d\ ' 


dt 2A, ' 


cZA A ._, dF 
~dT ~ ^\ t _, ' 


d\_, dF 

dt dA k _ 



Let us now consider the mean value of the function F relatively to the angular 



212 



G. W. HILL: EXTENSION OF DELAUNAY'S METHOD 



[April 



variables \ , \ 2 , • • • \ jt _ 1 . Then, since i'as a periodic function involves only- 
cosines of arguments, if [F~\ denote the mentioned mean value, we shall have 

(17) [F-\ = [i £J' W^\ 2 • • • dV, , 

where the first factor of the right member denotes an operation repeated k — 1 
times, once in reference to each of the variables X, , X 2 , • • • , \_ x • As F re- 
mains finite whatever values 6 , \ , X 2 , • • • , \ h may assume, it follows that [F~\ 
is finite whatever may be the value of 8 . Thus [F~\ is developable as a peri- 
odic function of involving only cosines ; and we may write 



(18) 



[F~\ = - B - A l cos 6 - A 2 cos 20 - A 3 cos 30 , 



where B , .4, , A 2 , A 3 , • • • are functions of the linear variables @ , A t , A 2 , 

•■•i A *-i- 

Let us now suppose that, in equations (16) , [F~\ is substituted iorF. They 

then become : 



(19) 



dt 

d\ 
'dt 



d'e" 



= o 



dA 



dt 



d6 
dt 

d\ 

dt 



d\,_ 
~dt 



d\FJ 

3JF-] 
dA,._, ■ 



A, , A 2 , • • • , A^, are therefore constants, and the two equations of the first line 
contain no other variables than ® and 6 , and thus form a distinct system by 
themselves and determine these two variables ; after which, by substitution of 



values, the remaining differential equations for \, , \ 2 , 
variables through quadratures. 



• i X,._, determine these 



As [_F ] involves only two variables O and , the two equations which 
begin (19) have the integral, C being an arbitrary constant, 



(20) 



[F] + C=0. 



This integral constitutes a relation between the two variables © and 6 ; and, 
if the latter are regarded as coordinates defining the position of a point in a 
plane, (20) is the equation of a plane curve. For this graphical exhibition of 
the connection between the two variables, we might adopt that in which they are 
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the polar coordinates of a point, © being the radius and 8 the angle. But, in 
some cases © may pass through zero. This difficulty may be obviated by adding 
to it a sufficiently large positive constant and thus it be rendered uniformly 
positive. This can be done provided it does not go to negative infinity. How- 
ever, all circumstances considered, it will probably be a better course to adopt a 
representation in rectangular coordinates, 8 being the abscissa and © the ordinate. 

If we derive from (20) an expression for 8 in terms of © and substitute it in 
the first equation of (19) and take the reciprocal of both members we shall have 
the time in terms of © by a quadrature, and, by the inversion of this, © as a 
function of t . On the other hand, if we derive from (20) an expression for © 
in terms of 8 , and substitute it in the second equation of (19) we shall have the 
time in terms of 8 by a quadrature, and, by the inversion of this, 8 as a func- 
tion of t . 

We proceed to note some of the properties of the curve whose equation is (20). 
In the first place it must be stated that if the differential equations of (19), 
which determine the variables © and 8 , compel the first of these to take on 
values rendering the right member of (18) a divergent series, we agree to set 
aside such cases as nugatory. Singularities of a certain kind are therefore ex- 
cluded. The curve cannot have a point d 'arret, for, at this point, we should 
have simultaneously d®/dt = and dd/dt = ; and in consequence all suc- 
ceeding derivatives of these variables would vanish. Thus, at this point © and 
8 would be invariable, which is impossible. It cannot have a multiple point, 
since, for given values of © and 8 , there is but one value of each of the quan- 
tities d®/dt and d8/dt . If the curve pass through a point, it must proceed 
thence until it returns to that point or goes on to infinity. In the latter case, 
taking a polar representation for the moment, it may either have two infinite 
branches, or may make an infinite number of turns about the pole, or, in other 
words be a spiral. But, since equation (20) involves only cosines of 8 without 
sines of the same, the curve must needs be symmetrically situated with respect 
to the axis from which 8 is measured. Hence, the last supposition must be re- 
jected ; that is, it cannot be a spiral, nor can it have more than one distinct turn 
about the pole. 

The curves graphically representing (20) may be divided into three classes. 
Here, for convenience, we adopt a rectangular representation. Let us suppose 
that an infinite number of values between the limits and it are substituted for 
8 in (20) ; the result will be an infinite number of equations for determining the 
corresponding values of © . Let one of these be satisfied by a real value of © . 
Then it may happen that all the remaining equations are satisfied by real values 
of this variable continuous among themselves and with the value first mentioned. 
The variable 8 can then move from — oo to + oo and there will always be a cor- 
responding real value for © . The first equation of (19) shows that © will beat 
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a maximum or minimum when = iir , i being a positive or negative integer. 
As, in the equation 

(21) - dt = A l sin + 2A 2 sin 20 + 3A S sin 80 + ... , 

the quantity A l is, in general, larger than A 2 , A z , ■ • • , it follows that © will 
have no other maximum or minimum values than those just mentioned. In ad- 
dition, if a maximum value occurs for = 2wr , then will a minimum value 
occur for = (2i + l)7r , and vice versa. If, in (20) we put, in succession 
= , = 7r , we shall have the two equations : 

\C-B= A l + A 2 + A i + ..., 
(22) 

[C-B=-A l + A 2 -A 3 + .... 

And if © be regarded as the unknown to be determined by them, it is plain that 
the maximum value of © will be a root of one of them and the minimum value 
a root of the other. Again © cannot be constant unless all the coefficients A 
vanish. It is quite evident that, in this case, the values of © and can be 
represented by the infinite periodic series : 

f © = © + ©, cos [0 O (< + c)] + © 2 cos 2 [0 O (< + c)] + • • ■ , 
(23) 

[ = e (t + c) + 0, sin [0 O (< + c)] + 2 sin 2[0„(< + c)] + • • .. 

These two equations are to be regarded as the integrals of the first and second 
differential equations of the group (19) ; c is one of the arbitrary constants in- 
troduced by the integration, the other may be supposed to be either the C of 
(20 ) or the © of the first of (23). But while C and c are conjugate to each 
other, this is not necessarily the case with the elements © and o (i + c) . The 
remaining coefficients of (23) , ©, , © 2 , • • • , O , l , 2 , • • • , are functions of C 
or © . On account of the form of the curve which represents (20) in this case 
it may be called the sinusoid case. 

We come now to consider the second case of the representation of (20) by a 
curve. Here, if we give to its range of values between and ir , we shall find 
that the equations determining the corresponding values of © have two real roots 
for an arc of values for 6 which either begins at or ends at it ; and, in the first 
case the arc terminates, or, in the second case, begins, at the same intermediate 
point. At this point the two real roots become equal, and, for the remainder 
of the semi-circumference, they are imaginary. Consequently, at this point, 
attains either a maximum or minimum value. Because the equation contains 
only cosines of multiples of , in the one case, the right line = 0, and, in the 
other, the right line = ir , divides the area embraced by the curve symmetri- 
cally. The maximum and minimum values of © are given by the roots of that 
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one of the two equations of (22) whieh has two real roots. In this case, 9 can 
not be represented by series like the second of (23), but, in general, we may 
give the integrals of the problem the form : 

f 6 cos = P + P x cos [6 a {t + c)] + P 2 cos 2[0 O (* + c)] + • • • , 

(24) { 

[© sin = Q x sin [0 o (t + c)] + Q t sin 2[0„(« + c)] + . . . , 

where O , _P , P x , P 2 , • • • , $i , Q 2 , • • • , are constant coefficients and func- 
tions of the C of (20) , while c , as before, is the other arbitrary constant. It 
will be perceived that, in the former case, the integral equations (23) can be 
given the form (24) if one chooses ; and Delaunay has always adopted it where 
the eccentricity e would appear as a divisor in the first form. At the two points, 
at which has attained its maximum or minimum value, we have d0/dt = , or 

,r»rN dB dA, . dA„ 

(25) -jK- + -T-=i cos + -^4 cos 20 + ■ • • = . 

v ' d® ^ d® d® ^ 

When dB/d® and dAJd® are quantities of the same order of magnitude, the 
second case is likely to occur. As the curve, which here represents the connec- 
tion of the variables ® and , is a closed one, this case may be called the ovaloid 
case. This kind of motion in the variables is, however, generally termed a 
libration. Observation has not yet shown that it occurs in the system of the 
eight major planets of the solar system, although it is possible it may exist for 
very large values of the integers j. . However, should this prove true, the in- 
fluence of this circumstance on the motion of the system would be quite in- 
significant. 

The third case in the graphical representation of (20) occurs when, in a 
certain range of values for , bisected by the value = or by the value 
= 7r , we find a real value for © , but this value tends towards positive or 
negative infinity as the limits are approached. Here there is one maximum 
and no minimum for ® or one minimum and no maximum. As in the previous 
cases, these values occur when = or = tr . As long as the instantaneous 
orbits of the planets composing the system are elliptic in their nature this case 
cannot present itself. And ® cannot go beyond a certain limit without some of 
the elements becoming imaginary. In order, therefore, to prevent the occur- 
rence of functions of complex variables, a modified system has to be adopted. 
But an illustration of this case can very easily be constructed. In order to 
escape the difficulty of divergence when |©| exceeds a certain limit, let us suppose 
that \F~\ is finite and does not run into an infinite series, and that all the quan- 
tities A i beyond A x vanish. Then the equation (20), being solved with reference 
to cos , gives 

cost/- ^ ~F(&)' 
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Let © be the value of © when 9 = ; in order to have the present case we 
ought to have 

a being less than unity. We may suppose that © is involved linearly in B and 
A 1 , so that, a , b , c , c? being constants, 

/(©) = a + 6© , i^(@) = c + d® . 
Then 

„ a — c 5 

All the conditions will be fulfilled if we put 

/(©)=3 + 4@, J P(©) = 2 + 5@; 

whence @ = 1 and a = |- . © is thus continuous while 6 is contained between 
the two values given by the equation cos 6 = | . At the limits © becomes in- 
finite. In a system of polar coordinates, if © is the radius and 6 the angle, the 
equation of the curve graphically exhibiting the connection of the variables © 
and 6 is : 

„ 3 + 4© n — 3 + 2 cos 6 

cos 6 = ^r^a , or © = 



2 + 5© ' ~ 4 - 5 cos 6 ' 

It is thus a quartic algebraic curve whose equation in rectangular coordinates is : 

[2sb - 4(x 2 + y 2 )f = (3 - 5x)\x 2 + y 2 ) , 
whose course resembles that of a hyperbola. The formula for the time is : 

Jc?@ /• d© 

If this be integrated between the limits © = 1 and © = © it will give the time 
required to describe the curve from the point 6 = to the point having the 
radius © . 

VI. 

Let us now suppose that by the integration of the system of differential 
equations (19), it is proposed to remove from F the periodic terms having the 
argument 6 , that is, those contained in [F~] . We confine ourselves to the first 
case as that will usually be the one which presents itself. The integrals of (19) 
will evidently have the form : 
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© =©„+©! cos [0 o (t + c)] + © 2 cos 2[0 O (« + c)] + ■ ■ • , 

A, = a constant , A 2 = a constant , • • • , A k _ " = a constant , 

= o {t + c) + 0, sin [0 o (t + c)] + 2 sin 2 [0 o (t + c)] + • • • , 

(26) X x = (X x ) + \?>(t + c)+ X« sin [0 O (« + c)] + Xf sin 2 [0 O (* + c)] + • • •, 

X 2 = (X 2 ) + X<°>(« + c) + X« sin [0 o (t + c)] + Xf sin 2[0 a (t + c)] + • • ., 

. V-i = OViH^+oH^ * in I?.(*+<0] +>& sin 2[^+c)] + • • -, 
where c , (Xj) , (X 2 ) , • • • , (X 4 _ 1 ) , A t , A 2 , • • • , A /t _ 1 , and one other involved in 
@ , © j , • • • , g , X ( j 0) , X ( j !) , • • • , "complete the number of 2k arbitrary constants. 
If then we take formulas of transformation as follows, replacing 

© by @ + ©j cos + © 2 cos 2 + ■ ■ ■ , 



(27) 



A, 



by A x , 



A k _, by A M , 

by + X sin + 2 sin 20 + • • • , 

Xj by Xj + Xf sin + Xf sin 20 + • • • , 

\^ by X,_ x + X« x sin + X« x sin 20 + ... , 

the terms of F involving as the argument will disappear from it, their mean 
values going over into the non-periodic term — B . 

We shall now have a new set of angular variables , X x , X 2 , • • • , \ h _ 1 , and 
it will be pertinent to inquire what are the linear variables conjugate to them, 
so that we may still have differential equations of the canonical form. Taking 
the two equations of (27) which define © and in terms of the new variables, let 
us suppose that the latter are © and 0' (an accent is applied to the latter for the 
sake of discrimination).. Then we have 

© = © + @ x cos 0' + @ 2 cos 20' + • • ■ , 

= 0' + y sin 0' + 2 sin 20' + • • • . 



(28) 



Then, by the substitution of these values in the two differential equations of (16) 
which determine © and , we get 



(29) 



' d® d® d® d0' 
d® W + W dt~ 


dF d® BF30' 
d® ~W + W Iff' 


B0 d® 30 d0' 
d® dt + W'dt~ 


dF d® dFd0' 
~d® d®~Wd®' 
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By multiplying these equations by the proper factors, and putting A for the 
functional determinant or Jacobian : 

d® B0 d® 30 

( 30 ) *~d9 t W~fflW ' 

we have 

d® A _dF ^__SF 

( dl ) A dt -dd" A dt ~ d® - 

But 

[d® d® ~\ r 1 

1 + ^ l cos $' + ^ 2 cos 20'+-. J 1 1 + ^cos 0' + 20 2 cos 20' + ■ ■ ■ J 

(32) 



+ 



[©! sin 0' + 2® 2 sin 20' + • • ] [^sin 0' + 1| sin 20' + . • •] . 



= 1 + 1 4, PA + 20 2 2 + 3^ 3 e 3 +...]. 



According to the theorem of Poisson, A is independent of t , or what in this case 
amounts to the same thing, of 0' . Hence, in computing its value, we have re- 
gard only to the absolute term. Thus 

A = x + 2 [^ + 2 ^e + ••• + ®i ^e o + 2e W o + ---J 

(33) 

_ i _l i _ . 
d® 

Then, if we adopt a new variable ©' in place of © such that 

(34) 0' = / A<?@ = O + i [^6, + 20 2 © 2 + 30 3 3 + •••], 

equations (31) will be transformed into : 

d®' _ dF d0' d_F 

which have the canonical form. As to the remaining linear variables A l , A 2 , 

• • , A t _, , which are identical with the former variables denoted by the same 

symbols, it is evident that they remain the conjugates of the new variables A t » 

\> " ' • ' \-i • 

VII. 

As it is somewhat difficult to discover the linear transformation of variables 
required to pass from the set 

L x , L 2 , • • • , L k 1 , A t , A 2 , • • • , A,_ 1 

t to the set 
h > h » • • • i h -I " » \ ' \ ' ' ' " ' ^t-i 
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suppose we adhere to the use of the first set. Then we have the differential 
equations 

d® _ d[F] d6 _ d[FJ 

~di ~ ~df~ ' dt~ d&~ ' 



(36) 



(37) 

Then 

(38) 



dL. _ . 3_[F] dl t 
dt ~ 3i ~W di 



d[Fj 

dL. ' 



= Wj^y (*"=1, 2, •••, k). 



A = (A)+J,@, 



« = 1,S, 



*), 



where the (L^) are arbitrary constants. By means of these equations the linear 
variables L. may be eliminated from the right members of all the differential 
equations (36) and (37), and thus be replaced by the single variable ® . Thus 
the mentioned right members become functions only of the two variables @ and 
; and, after © and are found in terms of t from the integration of (36), the 
L i will be found from (38), and the L from a quadrature operated on the second 
equations of (37). 

Then we can suppose that the integrals of (36) and (37) are written thus : 

0=6, cos [0 o (t + c)] + 2 cos 2 [0 o (t + c)] + • • • , 
A = (A) +ii°i cos [U* + c )] + Jfi, cos 2 W + c)] + • ■ • , 



(39) 



A = (A) +JA cos [6 (t + c)] +j k (® 2 ) cos 2 [0 o (t + c)] + • • • , 
= o (t + c) + 1 sin [0 o (t + c)] + 2 sin 2 [0 o (t + c)] + • • • , 
Z, = (Z,) + Zf (t + c) + Z<» sin [0 o (t + c)] + If sin 2[0 o (t + c)] + • • • , 

h = (h) + *T(t + c) + I? sin [0 o (t + c)] + Zf sin 2 [0 Q (« + c)] + • • • , 

where the (L^ , (L 2 ) , • • • , (L k ) , (Z,) , (Z 2 ) , • • • , (Z t ) are arbitrary constants. 
The latter set are not independent as they must satisfy the relation : 

As there is the additional arbitrary constant c , the number of independent con- 
stants is 2k as it should be. The coefficients O , l ( p , l ( }\ l ( P , • • • , which may 
be regarded as functions of the (L.) , satisfy the relations : 



W+JW + 
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The formulas of transformation are then : 

Keplace L x by L x + j i ® 1 cos 6 + j x ® 2 cos 20 + • • • , 



(40) 



replace L k by L k + j k ® x cos 6 + ,/ 4 2 cos 20 + 
replace ^ by l x + Z™ sin + Zf sin 20 + • • • , 



replace f„ by l k + Z« sin <? + If sin 20 + • • • . ' 

We have now to find what linear variables are conjugate to the new 
angular variables l x , l 2 , ■ ■ • l k . They are discovered immediately from a com- 
parison of (34) and (38) . As, from (34) , it appears that @ the arbitrary con- 
stant, which may be conceived as annexed to the series for © , is not the element 
conjugate to the angular element Q (t -f c) , but that the expression 

* [0 1 e 1 + 20,e, + 80 8 e a +...] 

must be added to it to produce the required conjugate, it is plain from (38) that, 
after the transformation, the new L i is no longer the exact conjugate of l i , but 
that we have for that element the value 

(41) L { + folOfr + 20 2 2 + 3^0 3 + • • •] (»=i, a, -, *) ■ 

VIII. 

In making one of Delaunay's transformations it is not absolutely neces- 
sary that we should employ the linear variables L i , which are the conjugates of 
the angular variables L , in the development of the various series needed ; we 
may use any others connected with the former by known relations. Then equa- 
tion (41) will inform us, at any stage of the transformations, what function the 
conjugates of the angular variables are of the used linear variables. Thus, in 
making his developments in the lunar theory, Delaunay has not used the ele- 
ments he calls L , G , IT, but has substituted for them others which he names 
a, e, j. 

Let us suppose that the new set of linear elements we determine to use are 
denoted by the symbols e x , e 2 , • • • , e k . Then, in order to form the equivalents 
for the dljdt, it will be necessary to know the values of the partial derivatives 
deJdL in terms of the e i . The number of these derivatives is Te 2 , and we shall 
have as many equations for determining them. Having the L. in terms of the 
e { , the general form of these equations will be : 

.„„. BL. Be, BL. Be„ 3L. Be, 

1 3 2 j k j 
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according as i and j are different or are the same. These equations divide 
themselves into h groups each containing h equations ; each group serving to 
determine the k partial derivatives of the e. with reference to one of the L. 
which we denote by L. . The functional determinant of each group of equa- 
tions is the same, being the Jacobian of the variables L l , L 2 , • • • , L h with 



reference to the variables e 1 , e 2 , 

shall have 

de, 

(43) 



Calling this determinant A , we 



dA 



dL. 

3 



(t) 



We shall much abbreviate the calculation of these partial derivatives by 
agreeing to compute dejdL. not directly but by a correction to be added to the 
value it had before the transformation was made. This correction will be de- 
noted by I)(deJdL) . Now the transformation being that made to remove from 
F all the periodic terms having the argument 

Jih+hh+ ••■+JA^ 

suppose that, in this transformation, we put 



(44) 
and also 



k= m®i + 20 2 ® 2 + H®a + • ••] , 



<-) -*>■£(£)♦£(£) 



+ 



+ 



dK 



&)• 



where (deJdL m ) denotes the value of the partial derivative before the transfor- 
mation ; then we shall have the group of linear equations : 



(46) 



dL 

~de. 



dL, 

de 



1 1, \dLj + de i JJ {dLj + 



+ 



+ •■■ + 



dL 

Be, 

dL, 

de,. 



\ m/ 



dL 

de 



Put now 

(47) 



x n \dLj + de 2 \dL m ) 



+ ■■■ + 



dL, 
de,. 



B 



m 



--'kK 

•Jh m 



+ ■■•+.?»- 



dA 



Then, evidently, 

(48) 



'(£) •(%) '<%) 



D 






P.K 
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Thus the corrections for each partial derivative, caused by a given transformation, 
are reduced to the product of three factors, one of which is a function of i , 
another a function of m , and the third independent of either. The elements 
e { , before the performance of any transformation, may be what we choose. 
Delaunay chose his a , e , 7 so that they were connected with L , G , H, 
the conjugates of the angular variables I , g , h , by the relations : 



L-y/pa, G=L^\-e\ H = G(l - 2y*) . 



Part II. 

Application of Delaunay's Method to the Minor Planet 
of the Hecuba Type. 

IX. 

Delaunay's lunar theory affords a plentiful assortment of the transform- 
ations just discussed, but their application in a case of planetary motion gives 
rise to more complex expressions. In the lunar theory it is possible to expand 
all coefficients in power series of all the parameters involved ; but, in a planetary 
theory where o , the ratio of the mean distances, is a considerable fraction, it is 
necessary to introduce the functions of a usually denoted by the symbol b['\ as 
also their derivatives with respect to a. It may therefore be profitable to give 
as simple an illustration as possible of these transformations where b[' } must be 
used. 

Let Jupiter be supposed to describe a circular orbit about the Sun, while a 
small planet, without mass, describes an orbit in the same plane. Let the radius 
of Jupiter's orbit be taken as the linear unit ; denote its longitude by e' + n't , 
and the masses severally of the Sun and Jupiter by ra and m . Let a , e , I , 
and g be the mean distance, eccentricity, mean anomaly, and longitude of the 
perihelion of the small planet, and r and v its radius and true anomaly. Put 



7 = g — e' — n't , L = <Sm a a , T = \^m a(l — c 2 ) . 
The function we have denoted by F will have the expression : 

TYt 

(49) F= 2" + »T + m'{[l-2rcos(u + 7) + r ! ]-*-rcos(u+7)}, 
where r and v are to be eliminated through the equations : 



r cos v = a (cos u — e) , r sin v = a VI — c 2 sin u , u — e sin u = I . 
The position of the small planet will be known when we know L , T , 1 , 7 . 



dL dF 


dl 


dF 


dt ~ dl ' 


di ~ 


~ dL' 


dT dF 


dy 


dF 


{ dt ~ dy ' 


dt ~~ 


~W 
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The differential equations for determining the latter are : 



(50) 



X. 

In order to give an illustration of the transformations named by Delaunay 
operations, let us select from the periodic development of F , which, from (49), 
plainly has the form ^A iit cos (il + i'y), all the terms having the argument 
= I -f 27 . These will be terms of long period in case the small planet has a 
mean motion nearly double that of Jupiter, which case has been extensively dis- 
cussed by astronomers, such a minor planet being called of the Hecuba type. 
Taking 8 as one of the angular elements, we see that we can adopt 7 as the 
other, and thus shall have 1= — 2y . In order to obtain \_F] from F we 
have the equation : 

m„ 1 r* m'dy 

(si) [*]-rf + ' T + ;I vi -«,.»(. + ,) + ,.■ 

remembering that r and v are now the same functions of —2y they were be- 
fore of I . The last term of F in (49) is here omitted as it contributes nothing 
to [F] . Put 

[1 - 2r cos (v + 7) + *■*]-* = l&v + &» cos (v + 7) 
(52) 

' +# 2 >cos2(t> + 7) + .-., 

where Z?^ is the same function of r that 6^° is of a . In order that this series 
may be convergent it is necessary that a (1 + e) < 1 . Let us put 

(53) A™ = - f"& 2 " cos i (2v -T)dl. 

Then we have 

" ~2a 



97h _____________ 

[F] = s-° + n'y/m^a (1 - e 2 ) + m! [f4<°> + _<»> cos 8 
(54) za 

+ A™ cos 20 + . . •] . 

The investigation will be facilitated if we now make a slight change in the 
dependent variables employed so that they have the following equivalents : 



© = s/m.a , T = s/m„a[2 - -s/1 - e 2 ] , 

(55) 

[ 0= 1+ 2g — 2e — 2n't, y = e' + n't—g. 
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Then the differential equations determining the formulas of transformation are : 

d[F~\ dd d[F] 



(56) 



~db 

dt 



dd 



= 



d6_ 

Hi 

dy 

~di 



d® 



Of these equations the integral of the third, T = a constant, furnishes the 
relation : 



(57) a = «[3 -V\ - e 2 ]- 2 = o[l - e 2 + 1 e" - \ e« + & e s ] , 

a being a constant. By means of this relation the variable a may be eliminated 
from [F ] which thus will contain but two variables, e and . The equations 
(56) have the canonical form, but we prefer to discard the variable T and to use e 
in its stead. Supposing then that [F] is made a function of the variables e 
and 6 , the differential equations for the latter are : 



(58) 



de 
~dt~~ 

de 

dt : 



1 V\ - e 2 [2 - VI - e 2 ] 2 d[F] 



s/m^a 



de 



1_ s/1 - e?[2 - s/1 - e 2 ] 2 d[F] 
e de, 



s/m a 



where the factor 



n/1 - e 2 [2 - v/1 - e 2 ] 2 



[l + xe 2 -^-^ 



Tf J e — • • • J • 



These equations form a group to be integrated by themselves. After this in- 
tegration is accomplished, 7 is derived through a quadrature of the equation : 



dj 
dt 



bt 



In this equation [i^] is a function of ® and T , but we have preferred to write 
it as a function of a and e ; thus : 



d\_F] _ d\F] da d[F] 3e 

~dT dV dT + ~de~ dT ' 



But 



© = 



\/m a 



2- Vl- 



T = s/m^a • 



v/1 - e 2 = 2 - s- . 
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Consequently 



da _ la Be 1 VI - e 2 [2 - <Sl - e 2 ] 



dr Vm.a 



Remembering that, with our adopted linear unit, a' = 1 , n' = s/m^ -\- m! , we 
have : 



1_ dy 1 T 0a B[F] x/l-e a [2- v/l-e 2 ] BJFJ 1 

n' dt ~ 17 ^7\ \_ Ba e Be J' 

But, adopting, as before, g for the longitude of the perihelion, this is more 
simply : 



1 tig 1 f BR x/l-e 2 [2- x/l-e 2 ] 3R1 



l =~ i [2a 

m \ L 



n' dt \ / m ' \ \ """ Ba e Be J 



J( 1 + S) 

2 



^ i 1 <$ 



3a ~ 2 — >/ 1 — e 2 w ' ^ 



(59) 1 



where 



<I( i O 



[2 _ s/ 1 _ e 2 ]" + 



2- v/1 



m' 1 /" 7 dfy 

R 



m c 



l p 



*/l — 2r cos (v — 7) 4- J" 2 



XI. 



In an application like the present, where the periodic developments of the 
various quantities are always tardily convergent, it is nearly impossible to give 
literal expressions for the coefficients. And, even if we consent to give to each 
coefficient its numerical value at once, the work of multiplying such periodic 
series together is very embarrassing, and the process easily leads to the commis- 
sion of errors. Hence we adopt the method of substituting for each quantity 
involved the special values of it at equal intervals in the motion of the indepen- 
dent variable through the semicircumferenee. With this method of treatment 
it is necessary to separate the cases of non-libration and libration. 

It is always an advantage in computation to have the members dealt with in - 
dependent of any linear and temporal units. To this end let us substitute for 
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the independent variable t , the variable t = e' + n't or the longitude of Jupiter ; 
also we put 

[J^ 7 ] m'. 



W-- 



rn n 



m„ 



The coefficients of the periodic development of W are then absolute numbers. 
The equations which, with (57), we shall use for the elaboration of the problem, 
are the three following equations : 



(60) 











w= 


a constant , 






ch 










e 






1 




V(l 


+ 


v\ a 




dd ~ 


Vl -6 


J [2- 


-VI- 


-e 2 ] 2 


de 


dg 


VC\ 


2 

-4- 


v\ c 


dR 

i da 




1 


as 




dr 


2- 


VI- 


>?dr 





V(l + v) a 



[2 _ Vl - e 2 ] 2 + 



W has the expression 



W= r [2 - ^1 - e 2 ] 2 + •(! + *)« 



2a 



1 _ e 2 + 2 >/l - e 2 
3 + e 2 



+ 






c?7 
"A 



(61) 



J. 

2a 



= 9-„[i + e 2 + K + ie 6 + A" 8 



* + rt 



1-clO 



+ ••■] 



+ x/(l + v) a [1 - e 2 + i^ - le 6 + - 6 y - T f*e M + • • •] + R. 

This equation contains as variables only e and 8 ; hence, since e should never be 
negative, the dependence of the two variables on each other may be shown graphi- 
cally by taking e as the radius and 9 as the angle in a system of polar coordi- 
nates. If we are given a pair of simultaneous values of e and 9 , it is obvious 
that by their aid we can determine the constant value of W '. Desiring to ascer. 
tain at what points on the axis the curve passes we make in (61) in succession 
6 = 0° and 9 = 180° and we get two equations of the forms : 



(62) 



D 

D: 



J/je + Jf 2 e 2 + Jf 3 e s + M^ + A 6 e 5 + 
Mfi + M/- — M B e s + ik 4 e* — Jf 6 e 5 + 



where JD may be regarded as the arbitrary constant and the M are constants, be- 
ing functions of a and v . These equations are transcendental in e and are such 
that the positive roots of the one are equivalent to the negative roots of the 
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other. If each has a positive real root continuous with the value of e which was 
used for the determination of the constant value of W, the variable generally 
moves through the whole range of real values. But, if the first equation has two 
positive real roots and the second none, will librate about the value 0=0°. 
But, if the second has two positive real roots and the first none, will librate 
about = 180° . It will be seen that when D = we have the limit separat- 
ing non-libration from libration. 

XI. 

Case I. — Non-libration. — Here, as 8 goes through the semicircumference, it 
can be employed as the independent variable. Then, in the first equation of (60), 
we assign to , in succession, a series of equidistant values covering the semicir- 
cumference. (Those used in our illustrative examples are 13 in number, viz, 
6 = 0° , = 15° , = 30° , • • • , = 180° .) This procedure furnishes us 
with a like number of equations for determining the corresponding values of e . 
Solving these by the tentative process we have these values of e , and can apply 
to them the procedure of mechanical quadratures. Thus is obtained a general 
expression for e as a periodic function of involving only cosines. 

As the next step these special values of e can be substituted in the right mem- 
ber of the second equation of (60). To the special values thus obtained for 
dr/d0 can be applied mechanical quadratures, and the resulting periodic series, in- 
volving only cosines of integral multiples of 6 , can be integrated with respect to 
this variable. This integral may be put in the form : 

(63) 6 o (t + c) = + /3 1 sin6'+ /3 2 sin20 + /3 3 sin3<? + .... 

Knowing O we are now in possession of the period of the inequalities we are en- 
deavoring to derive. The left member of this equation we shall designate as the 
time-argument, and, for brevity, denote it as f . In the next place we assign to 
fa series of equidistant values going from 0° to 180°, and, by a tentative pro- 
cess applied to (63), arrive at the corresponding values of . These correspond- 
ing values of can be substituted for in the expression of e as a periodic func- 
tion of , and thus we shall have the values of e whieh correspond to the equi- 
distant values of f . We can now readily derive the similar values of the two 
quantities e cos and e sin . To these we apply mechanical quadratures and 
thus obtain the periodic developments of these quantities in terms of f . 

As the last step in this work we can, through the last equation of (60), express 
dg/dr as a function of e and , and, by the substitution of the special values of 
the latter variables, obtain the special values of dg/dr which correspond to the 
equidistant values of f . To these apply mechanical quadratures and the periodic 
series for dgjdr is obtained. This being integrated we have the series for g , 
and the solution of the problem is completed. 
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XII. 

Case II. — Libration. — Here we are shut off from the use of 6 as an inde- 
pendent variable on account of its not going through the semicircumference. 
But this difficulty is surmounted by substituting for it another variable which 
does move continuously from — co to -f oo . In order to ascertain, in the case 
of libration, the limiting values of we have to solve the simultaneous equations : 

dW 
W = a constant , -=- = , 
de 

the unknowns being e and 6 . That is to say, a value of 6 must be found which 
will make the first equation have two equal roots for e . This can be done by 
a tentative process. If we assume 6 too large, generally, we shall not be able 
to discover real values for e from the first equation ; but, if 6 is taken too small, 
we get two values real but unequal for e . These two conditions must be brought 
as close as possible until we discover the point of passage from one to the other. 
In our illustrative example we escape the necessity of this tentative process by 
assuming as one of the two fundamental elements of the example not the D of 
(62) but the amount of libration. 

The amount of libration being thus either assumed or determined, let k denote 
the limiting value of sin 6 ; we then can put 

(64) sin = k sin yjr ; 

and the motion of yjr can be regarded as extending continuously from — oo to 
+ oo . Adopting the variable yjr for replacing 6 , the second equation of (60) 
takes the form : 

dr -p: — e k cos i|r 1 

(65) df = ^ {l + v > a v/T^7 2 [2 - •TT^ x/1-^sinV dW> 

de 

where the newly introduced radical must receive the sign of cos . We can 
now make yfr play the same role as 6 did in Case I, and there is need of no 
further explanations. 

XIII. 

We attend now to the integration of equations (60). The operation of De- 
launay's lunar theory which is numbered B3* has great affinity with that here 
detailed, and the two may be compared. He, it is true, has six variables to 
our four ; but, in comparing, his 7 should be made to vanish and his h then 
becomes indeterminate. 



* M6moires de 1'AcacUmie des Sciences, vol. XXVIII, p. 493. 
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The periodic development of the reciprocal of the distance between two planets 
as a function of the time has been given by Leverrier to terms of the seventh 
order inclusive, and those of the eighth order have afterwards been added by 
M. Bouquet.* We avail ourselves of this development and adopt the mode of 
Leverrier for noting the coefficients except in the portion which is a function of 
e alone. We put lA ( P = (1/j !) a 3 ' d 3 Up/da s , j = in the portion factored by 
cos 08 , j = 2 in the portion factored by cos , j = 4 in the portion factored by 
cos 26 and so on ; only the numerical factors are written since the A can easily be 
filled in as they always commence with A^ , and the lower index always increases 
by a unit in each step to the right. With Leverrier we put % for \e . This 
then is the development of a' /A , preserving only the terms involving the integral 
multiples of 6 as arguments : 



+ 



+ 



(66) 



+ 



+ 



+ 



5 5 s + | 4 +5-16-10]^ 



3 12 



-[4 + l]*+[l4 + j-— 6— 3]aM-[- 

, T271 203 .19 49 . m . 25 „ _H , 1 „ 
+ L^6-i44+4-8+ 20+ 2- 50 - 35 J^ J 0086 

[224-7 + l]^+[-f-f + 8 + 16+4]^+[ 1 f+f-49 
— 112+2+45 + 15]x 6 +[- 



8312 500212 280 , 1228 



45 



+ 0—236— 16 + 140+56l;t s j cos 20 

_[ 134 + f + 10+l]^+[ 4 |- 3 + 6 f+107-^-32-5], 5 

, T 64177 384789 12123 , 6249 , „, D 189 1nn „,"! ,l a 

+ [—w — 8cr-^cr+^r+ 348 -ir- 102 - 21 .F '} cos3e 



64177 384789 12123 , 6249 

8~ " 

275528 109972 



[ 2 -f- +f + 80+13 + l]^ + [ 



15 



15 



•1676 + 48 + 176 



, *. , ,.1 » , r 6259444 , 592976 , 234338 44756 10036 1432 
+ 54+6J;t 6 + | — ^ \- 



L 45 



15 



15 



+ 400 + 196 + 281 a: 8 j cos 40 



T 33797 , 50345 , 1795 , 245 , ,„ , ,1 . , T 5652235 , 9141589 . 210217 



+ 



36 ' 144 
48985 2020 775 



12 



24 



82—7 



>'} 



cos 50 



[ 18 f 16 + 7 -f 9 +4357+1024 + 174 + 19 + l]^+[- 



45378132 
35 



3695460 803616 150828 



+ 80 + 2156 + 720+116 + 81 x s | cos 69 



*Annales de l'Observatoire de Paris, vols. I, XIX. 
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, f 46064791 , 70738549 , 1880921 , 193921 . 4844 .469 „ . ) , „ a 
+ 1 252O- + ^2O— t -~W- + ^4™+-3- + ^- 22 - 1 h ,COS70 

f 552146674 . 213998824 2018552 . 926516 . 41380 7192 
+ ( 315 + 315 + 9 " l ~ 15 + 3 + 3 



+ 304+25 + 1 ^s cos 89. 

From this expression we must eliminate the variable a by means of its value 
in terms of e given by (57) . We put 



p= [2-v/l-e 2 ]- 2 -l. 

Let Leverrier's coefficient of cos^'0 be denoted thus : 

c Q A™ + vl«fl + c 2 A™ + • • ■ = £ cAf> , 

where c. is a function of e . Denoting the similar coefficient, after the variable 
a has been eliminated through (57) by ^fA^ j) , we evidently have 

A = (i + P )% + \ (i + py-'pc^ + *£^ (i + P y-* P v, + • • • • 

By means of this formula we obtain the following expression, in which a, the 
argument of the various quantities A { ?i\ is the constant a of (57) : 

•^ =g+ [0-l + l]z 2 + [0+ 0-8-9 + 3]^+[0+0+48+100 + 0-50+10]^ 

+ [0—4— 160— 600— 224 + 720 + 240— 245 + 35]** 

+ { -[4+l]*+[l4 + f + 2 -3]^+[_|115| 5 _?^ + 65 + 32-10]^ 

+ [f + !? + 6 f-*f-796+f + 190-35J4 cos, 

+ { [22+7 + 1] , 2+ [_^_^_ 56 + 4 + 4 ], 4+ [l_|_0 + 155 7+ 3 | 17 

+ 248-222-35 + 15], 6 +[- 8 f-^|3- 2 -?f 2 - 19 p 

+3680+2184— 496— 280 + 56~];t 8 | cos20 

+ { _ [ 134 +f + 10 + l],3+ [ 4 -f- 3 +!f- 5 + 559 + f -16-5^ 

+ [_^_^_^J3_26«9 +612 + 9 71 + 18 _ 21 J ;t , J cos3e 

+ { [f-Vf +80+13+l], 4 +[- 2 l| 2 - 8 -^0i- 1 -if -1068-48 

T6259444 , 1602224 , 1180994 , 630004 



+ 34 + 6]^+[- 



45 ^ 9 T '15 ' 15 
, 11372 7912 



800 + 28 + 28 1* 8 1 cos 40 
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, f T33797 , 50345 . 1795 . 245 , „ , ,1 , , f5652235 , 13594381 



36 ' 144 



, 469037 . 256585 , 4628 95 ro Jl , ) ,. 
+ | [^^4-^i?- 9 +4357 + 1024 + 174 + 194-l]^+[- 



45378432 



35 
5152104 1549888 473688 



19088—1704 + 240+88 + s\x s } cos 60 

, f 46064791 , 70738549 , 1880921 , 191863 . 4844 , 469 , M . , 1 , _. 
+ 1 25^ + "720-+~60-+^4--- i -"3- + -2- + 21 + 1 F TcOsW 

+ { "Sf 74 +™f 24 + 2 -gfg- 2 +^g+ 4 -fg + 7 f+304 + 25 + l}^ CO8 8g. 

In forming the value of dg/dr we need to know the derivative of the fore- 
going expression with respect to a . By noting the equation : 

dA m 

and changing our mode of noting the coefficients so that the number first given 
is the coefficient of A™ instead of Af J) , we have : 

a ~ = J+ [ — 1 + 0+3]^ 2 + [0— 16— 51— 24 + 15]**+ [0+96 + 444 + 400—250 

—240+70]^ 6 +[— 4— 328— 2280— 3296+2480 + 5760— 35— 1680+315]x 8 

+ { ~[5 + 2]^ + [^+ 49-3-12]^+[-^|y- 2 -P-23 + 388 + 110-60]^ 
+[^ + ^f + 3 -f _^_!f? + 1455+1085-280^ } cos • 
+ { [29+16+3]^+[ -lf- 6 -if -156 + 32+20]^+ [ 5 f* 

+ ^V 426l + 10 4_l 2 85-l 2 + 105]^+[-?f 6 - 2 ^ 

— 6 ill 6 ^_ ^. 8 + 29320 +10128— 5432— 1792+504] x s \ cos 29 

i f T 361 liio loq_l^1 . , r 282 77 , 16537 . 3795 . OQn inK T . 
+ J — I -g-+113+33+4j^+| — g- H j- -\ — ^ — 1-230— 105— 30J^ 5 

r 1528703 23961 21 189024 3 21543 11975 
+ L 80 40 40 2 + 2 

+ 3021—21—168]*' 1 cos 30 

+ { [350?+?| L 4 +279 + 56+ 5 J xt + [_?lW8_504O64_ 17608 _ 4464 

_ 7 0+ 24 o+42]^+[^ 5 ^+^^+ 1 -^+ 2 ^ 6 

+ ~^- _ 20624 — 5404 + 448 + 252 ] X s } cos iff 



(68) 
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+ f_ri85533 + ?^5 + 4|5 +554 + 85 + 6 1 5+ r 36203321 , 19222825 



144 ' 72 
1194659 , 293609 , 44855 



) , 293609 . 44855 .„„ ,._ _."| , 1 ,„ 
H g 1 g 633— 455— 56 #' > cos 50 



+ i [52023+^^ + 16143 + 4792 + 965 + 120 + 7] X s 



+[" 



71138952 73219472 6070728 22765 12 



35 35 



+ 2296 + 768 +-72~| ^ 8 1 cos 69 



-1 



587299425 , 93309601 . 4721157 , 230615 , 55475 . ,.„ . .„ . _> , „„ 
^040- + -360- + ^0- + -^- + -6~ +1539 + 161 +8^'cos79 

( 76 6145498 56 9296288 . 12872308 . 4533664 242860 
+ t ~315 h 315 + 15 + 15 + 3 + lb2m 

+ 2303 + 208 + 9 | X s cos 8». 

It is desirable to have the means of verifying these truncated developments of 
a' /A derived from the work of Levereier and M. Bouquet. In fact, by the 
application of the first of two following theorems, an error has been found in 
M. Bouquet's expression for (225) ; in the coefficient of JT 3 , — h should be 
substituted for h . The two theorems are the following : 

The coefficient of cosjO in the periodic development of a /A is the same as 
that of s> in the expansion of the expression 

5^(-x)<. + ^[i-x(. + i)](^)'^ 

in a power series with reference to s . 

The coefficient of cos j6 in the same development after a has been replaced by 
a [2 — v'l — e 2 ] 2 is the same as that of s 3 ' in the expansion of the expression 

| A*» [p - X (1 + p) (s + I)] [l - X (s + J)] (i^) ¥ «*W) 

in a power series with reference to s . 



In these expressions » stands for e/( 1 + v^l — e 2 ) . 

XIV. 

The two linear elements which determine all the coefficients in the periodic 
developments involved in this problem may be taken to be the constant a of 
(57) and the constant D of (62). It is proposed to elaborate two examples 
illustrating the subject in hand, one exhibiting non-libration, the other libration. 
In both we will assign to a such a value as makes log a = 9.8 . This value 
makes the period of revolution of the small planet nearly or exactly half that of 
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Jupiter. Whether we are to have a case of non-libration or Hbration will then 
depend on the value assigned to the second constant D . 

In the first place then we compute the values of such of the quantities A ( f as 
are needed in this investigation, correspondent to log a = 9.8 , by procedures 
which it is unnecessary to detail. The results are contained in the following table : 



o 

2 

4 



10 

12 
14 
16 





2 

4 

6 

8 

10 

12 

14 

16 



i = 
0.354 4041 774 
9.564 3962 993 
9.035 0709 047 
8.555 0516 205 
8.096 8549 86 
7.651 0634 5 
7.213 2942 8 
6.781 1495 
6.353 1544 



Values of log A { i ' ) for log a = 9.1 
i=l t=2 

9.845 4797 897 9.935 0116 655 

9.965 8367 1 0.002 7463 5 

9.694 9897 1 9.992 6030 6 

9.374 1611 5 9.842 7442 

9.031 7132 5 9.624 7440 8 

8.677 0326 9 9.367 4385 2 

8.314 4461 7 9.084 6837 2 

7.946 3186 8 8.784 1408 8 

7.574 0873 8.470 4851 



i = 
9.989 
0.007 
0.098 
0.099 
0.005 
9.849 
9.649 
9.419 
9.167 



= 3 

1230 

7852 

8192 

5969 

5654 

2650 

5418 

6593 

1668 



0.111 
0.121 
0.164 
0.231 
0.239 
0.178 
0.062 
9.905 
9.715 



4 

3716 

2342 

1565 

3986 

5776 

8047 



2677 
8220 



i = 
0.251 
0.257 
0.278 
0.329 
0.385 
0.398 
0.360 
0.272 
0.146 



= 5 

2555 

4382 

2372 

4253 

1308 

9502 

4400 

3635 

2268 



0.408 
0.412 
0.426 
0.453 
0.503 
0.552 
0.570 
0.545 
0.480 



= 6 

5603 

7086 

0622 

2441 

3035 

5313 

3163 

6897 

0088 



i=7 
0.576 748 
0.579 698 
0.588 970 
0.606 143 
0.637 049 
0.684 255 
0.729 206 
0.750 571 
0.721 959 



1 = 8 
0.753 37 
0.755 57 
0.762 42 
0.774 66 
0.794 30 
0.826 57 
0.870 86 
0.909 38 
0.935 78 



i = 9 
0.937 
0.938 5 
0.943 9 
0.956 8 
0.977 1 
1.004 5 
1.029 
1.045 5 
1.045 7 



Substituting these values in (67) we get 



A 



1.1307697497 + 
[ — 1.19571949 + 



0.04010033e 2 — 0.7367846e* + 1.17661e 6 +0.6155e 8 



(69) 



+ 



+ 

+ [1.70905245 

+ [ — 3.00445698 

+ [5.7966694 

+ [ — 11.800399 

-j- [24.86635 

— 52.40299e' cos 70 

+ 118.0918 e 8 cos 89 



3.1113902 e 2 — 
9.8883917 e 2 + 
27.190861 
- 71.99282 
+ 186.12652 
— 475.7506 



2.669146 
30.18579 
e 2 — 117.01214 
e 2 + 369.2943 
e 2 ]e 5 cos 50 
e 2 ]e 6 cos 6d 



1.90033e 6 ]e cos 6 
62.2O57e 6 ]e 2 cos20 



e«> 3 cos30 
e*]e 4 cos id 



"We adopt the mass of Jupiter so that v = 1/1047.355 . Then, in the ex- 
pression (61) of W, the portion which is independent of the interaction of 
Jupiter and the small planet, developed in powers of e 2 , becomes : 



1.58715 39467 862 
— 0.00226 07543 2 e 2 
+ 0.59490 01358 e* 
+ 0.09877 323 e s 

+ 0.13623 71 e 8 

+ 0.06804 8 e 11 



If we omit from the expansion of W its constant term, and call D the con- 
stant of the thus modified W, as iu (62), we have, as an integral of our problem, 
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Z>= — 0.00222246709e 2 + 0.5941966641e 4 + 0.09989661 e 6 + 0.1368248e 8 

+ [ — 0.00114165636 + 0.0029707121e 2 — 0.002548463e 4 — 0.00181441e 6 ]e cos 

+ [0.00163177054 — 0. 00944 12989e 2 + 0.02882097 e i — .593931 e 6 ]e 2 cob 20 

+ [ — 0.0028686138 + 0.025961456 e 2 — 0.11172156 e«]e 3 cos 30 

(70) + [0.0055345794 — 0.06873774 e 2 + 3525971 e*>* cos 40 

+ [ — 0.01126686 + 0.1777U0e 2 ] e 5 cos 50 

+ [0.02374204 — 0.4542400e 2 J e 6 cos 66* 

— 0.0500337e' cos 70 
+ 0.1127524e 8 cos 80 . 

By making = 0° in the preceding equation, we get, as the correspondent of 
the first equation of (62), the following : 



(71) 



— 0.00114 16563 6 e —0.00059 06875 5 e 2 
+ 0.00010 20983 5 e 3 + 0.59028 99445 e' 
+ 0.01214 61357 e 5 + 0.08372 1913 e 6 
+ 0.01414 140 e 7 + 0.08854 10 e 8 



= D. 



It will be seen by comparison of the coefficients of this equation that, unless 
e is very small, it will not do to regard the equation as approximately a quad- 
ratic in e ; for e = 0.1 the term in e i is ten times more important than the term 
in e 2 . The supposition that the mean motion of the small planet is nearly 
double that of Jupiter makes the coefficient of e 3 nearly vanish. In fact a very 
small change in the adopted value of a would make this coefficient . 

What sort of a curve we shall have exhibiting graphically the connection be- 
tween e and 8 will depend on the value assigned to D . To bring this out in a 
clear manner we compute the values of the left member of the preceding equa- 
tion for each 0.01 in the value of e between the limits ± 0.3 , and thus have 
the following table : 



e 


D 


e 


D 


e 


D 


e 


D 


0.30 


+ 0.0051 0216 


— 0.15 


+ 0.0004 5648 


+ 0.01 


— 0.0000 1148 


+ 0.16 


+ 0.0001 9224 


0.29 


44 8080 


0.14 


3 7472 


0.02 


2298 


0.17 


2 8623 


0.28 


39 2024 


0.13 


3 0675 


0.03 


3430 


0.18 


4 0095 


0.27 


34 1625 


0.12 


2 5066 


0.04 


4509 


0.19 


5 3896 


0.26 


29 6473 


0.11 


2 0477 


0.05 


5485 


0.20 


7 0298 


0.25 


25 6178 


0.10 


16715 


0.06 


6294 


0.21 


8 9590 


0.24 


22 0365 


0.09 


13659 


0.07 


6857 


0.22 


11 2071 


0.23 


18 8675 


0.08 


1 1166 


0.08 


7082 


0.23 


13 8067 


0.22 


16 0764 


0.07 


9115 


0.09 


6861 


0.24 


16 7912 


0.21 


13 6306 


0.06 


7400 


0.10 


6073 


0.25 


20 1960 


0.20 


11 4987 


0.05 


5928 


0.11 


4582 


0.26 


24 0580 


0.19 


9 6512 


0.04 


4622 


0.12 


— 0.0000 2236 


0.27 


28 4159 


0.18 


8 0599 


0.03 


3420 


0.13 


+ 0.0000 1129 


0.28 


33 3102 


0.17 


6 6982 


0.02 


2269 


0.14 


5695 


0.29 


38 7834 


0.16 


+ 0.0005 5411 


— 0.01 


+ 0.0000 1136 


+ 0.15 


+ 0.0001 1657 


+ 0.30 


+ 0.0044 8802 



As e ought always to be positive, in the first half of this table we may change 
the sign of e , provided we suppose that the corresponding value of D is re- 
garded as appertaining to the special value 180° for 8 , while, in the remainder 
of the table, this value corresponds to 8 = 0° . 



1900] 



TO THE GENERAL PROBLEM OE PLANETARY MOTION 



235 



From the course of the values of D in the table we see there is one minimum 
= — 0.00007 082, which occurs for e = 0.08 about; consequently, if D is 
chosen greater than this the equation (71) will have two real roots for e . If D 
is positive one of these roots will be negative ; changing the sign of the latter 
it will belong to the value 6 = 180° ; the positive root will belong to 6 = 0°. 
Thus, in this case, the motion of 6 is generally through the whole semicircum- 
ference, and hence is continuous from — co to + oo . But, if D is negative, 
both roots will be positive, and thus belong to 6 = 0° . In this case, therefore, 
departs from 0° and comes back to it without having reached 180°. This is 
called a libration ; we see that D = marks the dividing point between continu- 
ous and libratory motion for . The latter case also has the largest swing in the 
values of e , viz, from e = to about e = 0.127. Generally, the larger D is, the 
smaller will be the variation in e . Thus, if D — + 0.0045 , e will vary from 
0.29 to 0.30 . If there is libration e cannot exceed 0.127. These remarks, 
however, must be understood as applying only to the values holding for 0=0° 
and = 180°. Larger values for e may obtain for values of lying between 
0° and 180°. 

XV. 

For our illustrative example, in the case of a continuous motion for 6 , we 
assign to D the value + 0.0001 in (70) . All the coefficients of the various 
periodic series will now have determinate numerical values. The preceding 
table shows that, for this assumption, the eccentricity will have, when = 0°, 
the approximate value e = 0.1475 , and, when 6 = 180°, the approximate value 
0.0745 . In this case these are the limiting values, as e continuously 
diminishes while is passing from 0° to 180°. 

Attending now to the elaboration of our selected example, in (70) we give to 
, in succession, the values 15°, 30°, 45°, 60°, 75°, 90° and get as the right mem- 
ber of (70) , [the value 0=0° has already been considered in (70)] , 



15° 
30 


— 0.00110 27554e 
98 87034 


— 0.00080 93046e 2 
140 65774 


+ 0.00084 1072e 3 
257 2712 


+ 0.58878 755e« 
58670 872 


45 


80 72730 


222 24671 


412 9027 


58866 208 


60 


57 08282 


303 83568 


435 3970 


59615 002 


75 
90 


— 0.00029 54824 



363 56296 
— 0.00385 42466 


+ 0.00279 7293 



60514 036 
+ 0.60917 254 


8 
15° 
30 
45 
60 
75 
90 


+ 0.01297 98e 5 
+ 00755 03 

— 0.01219 27 

— 0.03286 91 

— 0.02990 01 




+ 0.09048 7e 6 
+ 0.12493 4 
+ 0.16863 4 
+ 0.14359 7 
+ 0.04056 8 
— 0.02140 4 


— 0.02180 7e' 

— 0.11214 3 

— 0.08332 4 
+ 0.17465 3 
+ 0.29851 4 




+ 0.20531e s 
+ 0.32869 

— 0.10302 

— 0.52039 
+ 0.30818 
+ 1.11581 



The coefficients in the second quadrant for are the same as in the first, but in 
reverse order, except, that for the odd powers of e the sign must be reversed. 
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Making the right members of these 13 equations equal to + 0.0001, we solve 
them with reference to e as the unknown, and substitute, in succession, the 
values thus obtained and the corresponding value of 6 in the right member of 
the second equation of (60) . The results obtained are the following : 




0° 


e 
0.14746 2372 


drjdv 
18.321384 


15 


0.14702 7366 


18.443969 


30 


0.14569 2650 


18.827993 


45 


0.14335 8947 


19.521595 


60 


0.13982 8670 


20.600475 


75 


0.13480 9503 


22.197678 


90 


0.12791 2165 


24.508178 


105 


0.11869 0066 


27.6150y7 


120 


0.10709 4867 


30.992715 


135 


0.09444 8152 


32.866681 


150 


0.08360 2031 


31.947061 


165 


0.07679 7488 


29.983003 


180 


0.07454 7767 


29.112462 



The mean of the numbers in the third column, attributing half weight to the 
first and last, is 25.101781 ; and this is the number of revolutions of Jupiter in 
the period of the inequalities we are investigating. If the sidereal revolution of 
Jupiter is put at 11.861980 Julian years, the latter period is 297.75681 such 
years. 

From the special values given in this table we can derive the two periodic 
series representing them. Integrating the latter, and, for brevity, putting % for 
& (t + c) , we get the following expressions : 



0.11918 
+ 0.03553 

— 0.00857 
+ 0.00123 
+ 0.00027 

— 00029 
+ 0.00012 

— 0.00001 

— 0.00000 
+ 0.00000 

— 0.00000 
+ 0.00000 
+ 0.00000 



891 

171 COS 

010 cos 
337 cos 
721 cos 
767 cos 
029 cos 
828 cos 
783 cos 
780 cos 
374 cos 100 
037 cos 116 
033 cos 120 



e 

■20 

w 

40 
50 
60 
70 
80 
90 



'+■ 



— 60173 

— 1643. 
+ 4612. 

— 2132. 
+ 544. 

6. 
87. 
52. 
16. 

0. 

3. 

1. 



+ 

+ 

+ 



40 sin 


e 


74 sin 


20 


77 sin 


30 


15 sin 


40 


20 sin 


50 


13 sin 


60 


51 sin 


70 


96 sin 


80 


90 sin 


90 


02 sin 


100 


89 sin 


110 


98 sin 


120 



The first of these is simply a transformation of the equation W = D by which 
e is expressed in terms of . From the second, by attributing to £ in succes- 
sion the 13 values 0°, 15°, 30°, • . • , 180°, using a tentative process, we can get 
the corresponding values of 6 . Thence by substitution in former results, the 
corresponding values of the four quantities e , e cos 6 , e sin 6 and dr/dO can be 
obtained. The results follow, the first column containing the argument : 
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? 
o° 

15 

30 

45 

60 

75 

90 

105 

120 

135 

150 

165 

180 



0° 

20 27 

40 25 

59 23 

76 57 

92 51 

107 3 

119 48 

131 35 

143 3 

154 48 

167 9 

180 



0.00 
58.05 
23.50 
15.52 
25.45 
50.00 
42.77 
15.87 
27.90 
46.20 
43.11 
46.91 

0.00 



0.14746 237 
0.14664 872 
0.14418 648 
0.13999 969 
0.13402 600 
0.12634 244 
0.11722 678 
0.10725 784 
0.09727 656 
0.08820 764 
0.08092 658 
0.07619 338 
0.07454 777 



+ 0.14746 237 
0.13739 211 
0.10976 569 
0.07129 163 

+ 0.03024 713 

— 0.00631 251 
0.03439 485 
0.05331 151 
0.06457 314 
0.07050 393 
0.07323 176 
0.07428 902 

— 0.07454 777 



e sin 6 

0.00000 000 
+ 0.05127 624 
0.09349 457 
0.12048 825 
0.13056 830 
0.12618 464 
0.11206 745 
0.09307 055 
0.07275 327 
0.05300 740 
0.03444 156 
+ 0.01692 847 
0.00000 000 



dTjae 

18.321384 
18.551258 
19.274084 
20.547277 
22.453723 
25.044275 
28.087856 
30.953203 
32.686384 
32.666482 
31.300251 
29.759410 
29.112462 



From the numbers in the fourth and fifth columns are derived the series : 



+ 0.00071 
+ 0.10563 
+ 0.03542 
+ 0.00539 
+ 0.00031 

— 0.00002 

— 0.00000 
+ 0.00000 
+ 0.00000 
+ 0.00000 

— 0.00000 

— 0.00000 
+ 0.00000 



143 

221 cos 
782 cos 
467 cos 
839 cos 
420 cos 
130 cos 
245 cos 
111 cos 
026 cos 
044 cos 
032 cos 
030 cos 



? 

2? 

3C 

4? 

5? 

6? 

7? 

8? 

9? 
10? 
11C 
12? 



+ 0.11737 
+ 0.03373 
+ 0.00528 
+ 0.00035 

— 0.00001 

— 0.00000 

— 0.00000 

— 0.00000 
+ 0.00000 
+ 0.00000 

— 0.00000 



247 sin 


? 


708 sin 


2? 


997 sin 


3? 


675 sin 


4? 


633 sin 


5? 


316 sin 


6? 


095 sin 


7? 


062 sin 


8? 


041 sin 


9C 


018 sin 


10? 


004 sin 


11? 



These forms for the integrals of our problem are to be preferred since they 
can also be used for the case of libration. 



XVI. 

To complete the solution the periodic series giving the position of the perihelion 
must be derived. Using logarithms^ instead of the actual coefficients, the first 
term of the right member of the third equation of (60) has the expression : 

. 2 g~= [6.9251786] + [7.126034V 2 — [8.00771]e J + [7.9475]e6+ [8.7115]e 8 

V (1 +v)a oa,- 

+ { — [7.9015032] + [8.284709> 2 + [7.66255]e 4 — [8.9051]e 6 } e cos B 

+ \ [8.3084076] — [9. 068286] e 2 + [9.55239]e* — [9.7559] e 6 }e 2 cos 20 

+ j— [8.707357] +[9.66685] e 2 — [0.2895] e*}e 3 cos30 

+ { [9.104403] —[0.20025] e 2 + [0.9097] e 4 }e*cos4e 

+ {—[9.50120] +[0.7008] e 2 }e 5 cos50 

+ | [9.89773] +[1.1811] e 2 }e"cos60 

— [O.2839]e'cos70 

+ [0. 6902] e 8 cos 80. 

The remaining terms of this expression for dgjdr can readily be derived from 
the values of e and dr/dd correspondent to the argument £ which have just been 
given. Calling these the second part of dgjdr , we have the following results :: 
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First Part 


Second Part 


dg\dr 


0° 


+ 0.00004 9796 


— 0.00579 1193 


— 0.00574 1397 


15 


0.00005 9848 


0.00572 9341 


0.00566 9493 


30 


0.00012 1787 


0.00554 3667 


0.00542 1880 


45 


0.00023 7477 


0.00509 7143 


0.00485 9666 


60 


0.00039 1638 


0.00420 4082 


0.00381 2444 


75 


0.00059 5969 


0.00270 5347 


— 0.00210 9378 


90 


0.00084 3915 


— 0.00024 2584 


+ 0.00060 1331 


105 


0.00109 9821 


+ 0.00326 3944 


0.00436 3765 


120 


0.00131 9503 


0.00773 7116 


0.00905 6619 


135 


0.00146 5221 


01284 0064 


0.01430 5285 


150 


0.00153 5401 


0.01789 1412 


0.01942 6813 


165 


0.00156 2180 


0.02178 6499 


0.02334 8679 


180 


+ 0.00156 9151 


+ 0.02328 0093 


+ 0.02484 9244 



The quantities in the last column furnish the periodic series for dg/dr . 
The absolute term shows that the mean motion of the perihelion of the minor 
planet is 0.00490 0079 times the mean motion of Jupiter. The integration of 
this series gives the expression for g . These two expressions follow ; (g) is the 
arbitrary constant added to complete the integral, and, in the second term, the 
unit of t is a Julian year. 

+ 00490 0079 

— 0.01441 6898 
+ 00444 6030 

— 0.00080 4772 
+ 0.00017 6760 

— 0.00006 7470 
+ 0.00003 1972 

— 0.00000 2614 
+ 0.00000 0139 

— 0.00000 3981 

— 0.00000 1706 
+ 0.00000 0414 
+ 0.00000 0648 

It is of interest to know the mean motion of the small planet which is not 
obvious at the beginning of the solution. We have the equation : 



dg 
dr" 



cos 


? 


COS 


2? 


COS 


3? 


COS 


4? 


COS 


5? 


COS 


6? 


COS 


7? 


COS 


8? 


COS 


9? 


COS 


10? 


COS 


11C 


COS 


12? , 



g) + 535.3662 


t 


— 74645.14 


sin ? 


+ 11509 92 


sin 2? 


— 1388.93 


sin 3? 


+ 228.80 


sin 4? 


— 69.87 


sin 5? 


+ 27.59 


sin 6? 


— 1.93 


sin 7? 


+ 0.09 


sin 8? 


— 2.27 


sin 9? 


— 0.88 


sin 10? 


+ 0.19 


sin 11? 


+ 0.28 


sin 12? 



d (I + g) dd dg 



d-r 



dr 



+ 2. 



Substituting in the right member the mean motions of 6 and g, its value is 
found to be 2.03493 7731 ; then, if for Jupiter we have /*' = 299."12838, for 
the small planet /* = 608.' / 70762. 



XVII. 

Illustration in the Case of Libration. — In the example we have chosen 
to illustrate the theory, libration, when it exists, is always about the value 
& = 0°. In addition to the value log a = 9.8 let us assume that the D of (70) 
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is to be so chosen that the half-swing of 6 may be 50°. Making, therefore, 
8 = 50° in (TO), we get the first of the following equations in e, and the second 
by taking the derivative of the first with respect to e : 



— 0.00073 38426 


e 


— 0.00250 58226 


e* 


+ 0.00439 3829 


e 3 


+ 0.58063 532 


e 4 


— 0.02026 791 


e 5 


+ 0.17135 53 


e<> 


— 0.01446 68 


e> 


— 0.32494 7 


e 8 



= £>, 



— 0.00073 38426 

— 0.00501 16452 e 
+ 0.01318 1487 e- 
+ 2.36254 128 e 

— 0.10133 955 e 
+ 1.02813 18 e 1 

— 0.10126 76 e 1 

— 2.59957 6 e 



= 0. 



Both of these equations should be satisfied when is at the limit of its swing, 
viz., when = ± 50°. The root of the second equation which is applicable to 
our purpose is e + 0.07606 124 , and this value substituted in the first gives 
D = — 0.00048 63102 , which is the value of D which brings about a libration 
of 50°. 

From the equation sin 9 = sin 50° sin yfr we obtain the following correspond- 
ing values : 












0.00 


15 


11 


26 


8.27 


30 


22 


31 


15.64 


45 


32 


47 


51.90 


60 


44 


33 


38.75 


75 


47 


43 


35.34 


90 


50 





0.00 



By means of these values we determine the form of (70) corresponding to the seven 
values of i^ . The coefficients are given in the following table (the small figures 
at the top of the columns denote the order of the final decimal) : 



i> 


e 


e 2 


e 3 


e* e 5 


e 6 


e 7 e s 


0° 


—114 16564 10 


— 59 0688 9 + 10 210 s +59028 994' +1214 61' + 8372 2 6 + 1414 5 + 8854 5 


15 


111 89923 


71 8977 


54 217 


58936 012 1284 01 


8713 9 


— 640 15961 


30 


105 45925 


106 9471 


164 931 


58751 743 1188 52 


10353 7 


6636 30349 


45 


95 96627 


154 8255 


291 587 


58665 083 + 487 32 


13434 5 


12360 30000 


60 


85 42473 


202 7039 


385 516 


58769 012 — 671 35 


16174 8 


11285 + 4731 


75 


76 79586 


237 7533 


429 468 


58965 923 1663 22 


17134 1 


— 4928 —22703 


90 


— 73 38426 


-250 5823 


+439 383 


+59063 532 —2026 79 


+17135 5 


+ 1447 —32495 



The expressions in this table constitute the left members of 7 equations of the 
8th degree ; they must be equated to the same quantity D = — 0.00004 863102. 
The two smallest real roots of each should be derived (they are those suited to 
our purpose). The connection of these roots with the variable ijr is settled in 
following way : the larger of the two roots is made to correspond to the value 
of ijr standing as the argument in the table, while the smaller is assigned to the 
value 180° — y}r ; the two roots being equal for yjr = 90°, the common value is 
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assigned to that value of ty . This arrangement is made in order that ty and r 
may augment together. These values of e together with the corresponding 
values of are, in succession, substituted in (65) ; thus we have the values of 
dr/d-fr corresponding to equidistant values of -yjr . These results are contained 

in the following table : 

f e drjdf 

0° 0.10846 187 37.23986 

15 0.10765 795 37.70988 

30 0.10518 818 39.14539 

45 0.10088 191 41.57020 

60 0.09452 189 44.76280 

75 0.08605 349 47.65995 

90 0.07606 124 48.00164 

105 0.06601 592 44.40759 

120 0.05744 757 38.68685 

135 0.05102 653 33.34281 

150 0.04671 683 29.47279 

165 0.04426 767 27.22320 

180 0.04347 566 26.49213 

It should be noted that, in the computation of the third column, the factor 
cos i]r/(d W/de) takes on the indeterminate form 0/0 ; employing the usual me- 
thod of treating vanishing fractions, this factor equals — 1/ \{d 2 Wjde 2 ) (de/d-*]r)~] . 
If here we should use the equation W=D to determine de/dty, the result 
would again be indeterminate. But this difficulty is avoided by employing the 
value of e as a periodic function of ty given by the quantities of the second 
column. Thus if 

e = « + «j cos f + «j cos 2-yjr + a 3 cos 3i/r + • • • , 

then 

de 

~jY = — a i sin ty — 2a 2 sin 2i|r — 3a 3 sin 3->^ — . • • , 

and, for the special value ty — 90° , this becomes : 

de 

djr = ~ a i + 3a s ~ K + " a 7 • 

For this special value of -^ it is found that 

d 2 W dp 

e ~ = + 0.00288 92836 , j— = - 0.03939 373 . 

The mean of the numbers in the last column of the table, allowing half weight 
to the first and last, is 38.65409, which is the number of revolutions of Jupiter 
contained in the period of libration ; thus this period is 458.5144 Julian years. 
From the special values of e and drfd-^r given in the table we derive the periodic 
series representing them. The latter can be integrated, and, as before, we put 
% for d (t + c) . Thus we get the following expressions : 
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0.0759 
+ 0.0338 

— 0.0000 

— 0.0015 
+ 0.0000 
+ 0.0001 

— 0.0000 

— 0.0000 
+ 0.0000 
+ 0.0000 

— 0.0000 

— 0.0000 
+ 0.0000 



8399 
8957 cos 
4144 cos 



2* 



2988 cos 3* 



4* 
5* 



3031 cos 
5016 cos 
0467 cos 6* 
1929 cos 7* 
0062 cos 8* 
0281 cos 9* 
0012 cos 10* 
0026 cos 11* 
0008 cos 12* 



c=*+ 



+ 29862. 37 sin * 




— 20922.46 sin 2* 




— 416.31 sin 3* 




+ 1650.00 sin 4* 




+ 13.01 sin 5* 




— 193.31 sin 6* 




— 0.37 sin 7* 




+ 27.39 sin 8* 




+ 0.01 sin 9* 




— 4.65 sin 10* 




— 0.02 sin 11* 




+ 0.81 sin 12* 




equation W = D , by which 



The first of these is simply a transformation of the equation W ■■ 
e is expressed in terms of the auxiliary variable yjr . 

Attributing to f , in succession, the 13 values 0° , 15° , 30° , • • • , 180° , 
by a tentative process we can get the corresponding values of ^ , as also, by 
substitution those of e , e cos 6 and e sin 6 . These results follow : 



s 


* 


e 


e cos 


e Bin o 


0° 


0° o' 0.00 


0.1084 6187 


+ 0.1084 6187 


0.0000 0000 


15 


15 30 0.64 


0.1068 0283 


0.1045 4084 


+ 0.0218 6455 


30 


30 35 58.93 


0.1022 4969 


0.0941 5541 


0.0398 7176 


45 


44 58 46.75 


0.0956 9588 


0.0804 5265 


0.0518 1767 


60 


58 28 33.58 


0.0879 7453 


0.0666 2875 


0.0574 4671 


75 


71 8 20.78 


0.0797 1930 


0.0549 1390 


0.0577 8952 


90 


83 15 23.39 


0.0714 4655 


0.0463 7246 


0.0543 5260 


105 


95 20 46.24 


0.0636 2489 


0.0411 4864 


0.0485 2749 


120 


108 7 15.59 


0.0566 9409 


0.0388 6515 


0.0412 7614 


135 


122 25 27.98 


0.0510 7994 


0.0389 6451 


0.0330 2919 


150 


139 5 9.66 


0.0468 5171 


0.0405 2863 


0.0235 0559 


165 


158 31 3.57 


0.0443 2198 


0.0425 4216 


+ 0.0124 3393 


180 


180 0.00 


0.0434 7566 


+ 0.0434 7566 


0.0000 0000 



From the data of the fourth and fifth columns result the periodic series : 





+ 0.0604 2349 




+ 0.0309 3374 cos s 




+ 0.0147 2940 cos 2s" 




+ 0.0015 8319 cos 3s" 




+ 0.0007 3464 cos 4s" 




— 0.0000 1013 cos 5s 


€ COS # = 


+ 0.0000 6769 cos 6s 




— 0.0000 1229 cos 7s" 




+ 0.0000 1611 cos 8s" 




— 0.0000 0669 cos 9," 




+ 0.0000 0106 cos 10s 




+ 0.0000 0528 cos lis" 




— 0.0000 0363 cos 12C 



! sin ft - 



+ 0.0571 8416 sin 


s" 


+ 0.0093 8541 sin 


2; 


+ 0.0030 1596 sin 


3; 


+ 0.0000 5256 sin 


4s" 


+ 0.0001 9928 sin 


5; 


— 0.0000 1597 sin 


o; 


+ 0.0000 1768 sin 


7s 


— 0.0000 0390 sin 


8s 


+ 0.0000 0771 sin 


9: 


— 0000 0714 sin 


10s 


+ 0.0000 0492 sin 


11C 
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XVIII. 

In computing, for this case, the values of dg/dr by the third equation of 
(60) we make a like division into two parts as in the former case. Substi- 
tuting the values of e and 6 which correspond to the values 0° , 15° , 30° , • • • , 
180° of f we get the following special values: 



? 


First Part 


Second Part 


dgjdr 


0° 


+ 0.00019 


583 


— 0.05008 


981 


— 0.04989 398 


15 


20 


725 


4843 


252 


4822 527 


30 


24 


656 


4385 


892 


4361 236 


45 


30 


013 


3722 


995 


3692 982 


60 


36 


434 


2938 


834 


2902 400 


75 


42 


685 


2098 


623 


2055 938 


90 


47 


979 


1247 


567 


1199 588 


105 


51 


817 


— 0.00410 


869 


— 0.00359 052 


120 


54 


013 


4- 0.00397 


237 


+ 0.00451 250 


135 


54 


781 


1211 


642 


1266 423 


150 


54 


497 


1991 


948 


2046 445 


165 


53 


746 


2629 


903 


2683 649 


180 


+ 0.00053 


360 


+ 0.02891 


583 


+ 0.02944 943 



From the quantities in the last column we obtain the periodic series for dg/dr , 
and thence by integration the expression for g ; these results follow : 



dy_ 

dT 



— 0.0116 3174 






— 0.0373 2443 


cos 


s 


+ 0.0007 9929 


COS 


2? 


— 0.0020 0735 


COS 


3? 


+ 0.0005 0972 


COS 


4? 


— 0.0002 9143 


COS 


5? 


+ 0.0000 8500 


COS 


6C 


— 0.0000 4117 


COS 


7? 


+ 0.0000 1080 


COS 


8? 


— 0.0000 0438 


COS 


9? 


+ 0.0000 0251 


COS 


10? 


— 0.0000 0293 


COS 


llC 


+ 0.0000 0214 


COS 


12? 





(*)- 


- 1270.844 l 




— 297586.8 


sin ? 




+ 


3186.4 


sin 2? 




— 


5334.9 


sin 3? 




+ 


1016.0 


sin 4? 




— 


464.7 


sin 5? 


9 = 


+ 


113.0 


sin 6? 




— 


46.9 


sin 7? 




+ 


10.8 


sin 8? 




— 


3.9 


sin 9? 




+ 


2.0 


sin 10? 




— 


2.1 


sin 11? 




1 + 


1.4 


sin 12? 



The unit of t in the second expression is a Julian year. 

By using the same formula as in the former case we find that the mean fi of 
the small planet in this case has the value 609".47474 . 

XIX. 

In attempting to apply the preceding method to the case where D = , we 
should find that dg/dr became infinite at the point where 6 = 0° ov d = 180°, 
and, when D is quite small, we should have to deal with inconveniently large 
numbers. This difficulty is surmounted by computing the differentials of the 
two quantities e cos g and e sin g in place of that of g . 



